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Abstract 

Casini et al raise the issue that the entanglement entropy in gauge theories is ambiguous because 
its definition depends on the choice of the boundary between two regions.; even a small change in the 
boundary could annihilate the otherwise finite topological entanglement entropy between two regions. In 
this article, we first show that the topological entanglement entropy in the Kitaev modeb which is not 
a true gauge theory, is free of ambiguity. Then, we give a physical interpretation, from the perspectives 
of what can be measured in an experiement, to the purported ambiguity of true gauge theories, where 
the topological entanglement arises as redundancy in counting the degrees of freedom along the boundary 
separating two regions. We generalize these discussions to non-Abelian gauge theories. 
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I. INTRODUCTION 


Casini et al explores the issue of defining entanglement entropy for gauge theories. Due to gauge 
redundancy, it is subtle to decide how to split the Hilbert space into products of subspaces each 
containing states within a given spatial region.^ To make the question precise, Casini chooses to 
dehne a Hilbert space by the operator algebra that acts on the space. Having chosen the algebra 
Av defining the Hilbert subspace of region V, its commutant My = Ay is thus the operator algebra 
that defines the Hilbert subspace of the complement region V. It is generally possible that Ay 
and Ay share a common center C, in which case the Hilbert space cannot be split into a product 
of two subspaces acted on by Ay and Ay separately. When such a common center is non-trivial, 
the computation of the entanglement entropy between V and V has to be properly dehned. It is 
argued that the most natural way is to express the density matrix p in the basis such that the 
common center C is diagonalized. One then throw away all off-diagonal components in p that maps 
states with different eigenvalues with respect to C, and each of the diagonal blocks characterized 
by given eigenvalue of C would now have a clear product structure of states acted on by Ay and 
Ay, from which we can obtain blocks of reduced density matrix by tracing out states in Ay. The 
final entanglement entropy would thus take the form as a sum of two parts, the classical part and 
the quantum part, namely, 

Se =-^Pilnp*-tr(pylnpy), (1) 

i 

where pi is the overall coefficient of a given diagonal block in the density matrix p. 

It was questioned how the choice of Ay affects the value of Se, particularly in the case where the 
ground state concerned is a topological phase characterized by long-range entanglement, and with 
a distinct topological entanglement entropy in the case of 2-1-1 dimensional topological ordered 
phases.This question is explored specifically in the context of 2-1-1 dimensional lattice gauge 
theories. In the specific examples of Kitaev models corresponding to Abelian groups, it is shown 
that one can readily find choices of Ay for which the entanglement entropy dehned above is 
identically zero, and that subsequently the topological entanglement entropy is also zero. We 
would like to elucidate the subtleties in this problem. 


II. ENTANGLEMENT ENTROPY IN THE KITAEV MODEL 

Let us verify that the Kitaev models have the usual topological entanglement entropy, quite 
independent of the choice of Ay. 
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A. Electric vs magnetic center in the Kitaev model 


The Kitaev model defined in Ref.l is not truly a gauge theory, in the sense that there is not 
any gauge redundancy in the Hilbert space, except that the ground state is chosen to be within 
the gauge invariant subspace, satisfying 


K\'^) = IV'), (2) 

where A^. are the star operator defined as in Ref.l. 

K =n • (3) 

l&Vi 

where L^i acts on the link variable on the link I and multiply the link variable 
element g, and I G Vi refers to all the links meeting at the vertex Uj. Whether it is 
multiplication depends on whether an oriented link I points toward or away from 
The convention chosen in Ref.l is such that 

Lf|h)in = |5'/i)in, Tf|/i)out = \hg~^)ont (4) 

and the state on a link |h)in,out is defined where h E G. 

The idea of the electric center and the magnetic centers are basically two natural choices of Ay 
in a lattice gauge theory with gauge group G. Consider a set of vertices V = {vi}. The electric 
centre corresponds to choosing a rough boundary for region V, namely that one includes all the 
operators Lf generating group multiplication on the link variables I that connects to at least one 
vertex v in region V. 


by the group 
a left or right 
the vertex Vi. 



FIG. 1. (Color online) Electric (rough) boundary on the left and magnetic (smooth) boundary on the 
right. The boundaries are in dashed lines, which are part of region V. 
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B. Naive entanglement entropy with the magnetic center 


Before we move on to computing entanglement entropy based on the choice of the operator 
algebra, let us consider computing the entanglement entropy in the most naive way, where we 
separate the links into 2 sets, without cutting links as being usually considered to preserve the 
original lattice structure. Consider in particular a convex region surrounded by a rough boundary. 
We recall that the ground state of the Kitaev model is given by 

(5) 

Vi g 

where |0) is the ‘empty’ state in which the degree of freedom on each link is set to the identity 
element of the gauge group G, and the factor J\f in the front is to ensure correct normalization. 
We have 

AA=|G|"^/l 

We can then rewrite the above ground state as follows 

IV.) = ® ife}) ® p) 

i,j,k 

where {q} denotes the set of conhgurations of the links in region A, the complement of region A, 
{aA is the set of conhgurations inside region A except for links along the rough boundary, and 
{hi} conhgurations inside the rough boundary. Note that i, j, k can be composite indices, as to be 
seen soon. Coefficients = 1 if a conhguration is allowed, and zero otherwise. 



FIG. 2. (Color online) An infinite lattice divided into two regions: A, inside the red (thicker line) loop, 
and A, ontside of the loop, the blue (dashed) lines comprising the rough boundary also belong to A. The 
wiggly lines simply represent the parts of the infinite lattice that are neglected. 
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Consider the lattice in Fig. 2, which is divided into two regions A, inside the red (thicker) loop 
and including the rough boundary (blue dashed) links, and A. Suppose there are exactly L links 
comprising the rough boundary, and the numbers of sites outside and inside the red (thicker) loop 
are respectively m and n, and let M = and N = We can rearrange the ground 

state as follows (and omit the curly brackets when necessary to avoid clutter): 


= l{ci}> 0 ^ (|6?.> 0 l{a,}) + \K,) 0 Kaj}) + ■ ■ ■ |6?„> 0 !{<■«})) 

- 9 

+I{C2}> 0 (|(.?.> 0 |{ai}) + m 0 |{02}> + ... |6f„> 0 |{a„})) 

. 9 

+1(^3}) 0 (|6?1> 0 |{a,}) + m 0 |{a,}> + ... |6|„) 0 |{a«})) 

. 9 

+ |{C|Gh-i}) ® + ■ ■ ■ ® 


+ 


L 9 


( 8 ) 


where each |{cj}) does not correspond to a single conhguration but rather the sum of \G\'^~^ 
different conhgurations in A, and each |{ai}) corresponds to the sum of different conhg- 

urations inside region A. It is clear that there are precisely terms in the above equation, 

in accordance with Eq. (5). 

The reason it takes the above form is that the ground state contains all linear combination of 
possible ways of acting on |0), leading to \G\'^^'^~^ different conhgurations of G elements on 
the lattice links. We can group all possible conhgurations of the links in A by the conhgurations 
on the links on the rough boundary. The L links on the rough boundary give rise to exactly \G\^~^ 
boundary conhgurations, each of which is associated with precisely conhgurations in A 

because there are m — L sites in A. Therefore, each |{cj}) is the sum of \G\^~^ conhgurations in 
A that share the same conhguration of the rough boundary. This way, the conhgurations inside 
the red loop can be considered separately. Now for each |{cj}), the group elements on all the links 
on the rough boundary are hxed at once; hence, regardless of the detailed group values on these 
links, there would be exactly |G| possible simultaneous choices of these values, resulting in the 
\bij) in Eq. (8). Consequently, for each rough boundary conhguration, namely for each |6y), 
we can similarly group the conhgurations in A excluding the rough boundary by the actions 
on the sites v E A that are directly at the end of the links on the rough boundary, and there are 
just L — 4 such sites. Thus, we have N groups of |{ai<j<Af}) conhgurations, each of which is the 
sum of conhgurations in A excluding the sites right at the rough boundary, the number 

of which is exactly n — L + A. 

If we now obtain the reduced density matrix, in which we trace out |ci), the vectors inside 
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the big bracket are clearly orthogonal eigenstates of the reduced density matrix. Since there are 
exactly such eigenstates, and that each orthogonal eigenvector can be normalized and the 

eigenvalues sum to 1, the von-Neumann entropy is 

SE = ln\G\^-^ =Lln\G\-ln\G\. (9) 

We thus recover the topological entanglement entropy 7 = In |G| which is the quantum dimension 
of the Kitaev model. 

In the Kitaev model, had we chosen to have smooth boundary instead of rough boundaries, 
the same arguments can be run through using the projector operator Bg. Since there is exact 
electric-magnetic duality, the result is independent of such choices of boundaries. 

When there are concave boundaries, then the number of vertices outside of A that controls 
the rough boundary would be less than the number L of boundary links on the rough boundary. 
However, that number depends on the precise form of the boundary and the number of vertices 
around the rough boundary outside of A scales linearly with L. Therefore, in the thermodynamic 
limit, we have a stable answer for the topological entanglement entropy that is truly L independent. 


III. ENTANGLEMENT ENTROPY IN “TRUE” GAUGE THEORIES VS SPT 

As noted in the previous subsection, we see that the topological entanglement entropy does not 
depend on the electric/magnetic boundaries for the Kitaev model. We would like to understand 
here why the calculation in Ref. 2 has led to vanishing topological entanglement entropy for mag¬ 
netic centers and in fact the same is true for a maximal tree structure in a true gauge theory. A 
key characteristic in the representation of the Hilbert space of a gauge theory is the redundancy 
in the labeling of the states: two conhgurations related by a gauge transformation generated by 
(3) are one and the same. Instead of claiming that the ground states are invariant under (3) by 
summing over all the states generated by (3) acted on different number of vertices, one requires 
that each orbit of the gauge transformations corresponds to one single state. 

Let us take the most extreme example in which we consider a single spin 1/2 system, in a 2d 
Hilbert space. 

^ = 1;)}. (10) 

Now suppose we just add a redundant label to the above states 

^ = {1 t,t), I(11) 
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which clearly is completely trivial. Consider the state 


|i'> = 4it,t) + 4ir4.). (12) 

Let us take as Pauli matrices acting on the “first” spin. While there is not really a second 
spin degree of freedom, we can always make the claim that that acts on the “second spin” takes 
eigenvalues satisfying 

a,' X = 1, (13) 

which is the direct analogue of gauge conditions in a gauge theory, although in a much more 
transparent form here. Now using basis diagonalizing the density matrix is 


P= l^)(^l 


1 

2 



(14) 


which is a two by two matrix, unlike the density matrix of a true 2 spin systems, which would have 
been a 4 by 4 matrix. In this case, is there a sense in which there is any entanglement? Clearly, 
there is not. However, entanglement entropy could also be viewed as the ignorance accompanying 
the dumping of some operators from the set of observables which we measure, i.e. the inability to 
determine the precise form of the density matrix if we only have measurements of a subset of all 
the observables. So one could talk about the density matrix, mostly likely mixed, with maximal 
entropy that reproduces the measurements of the remaining observable. So for instance, if we are 
not allowed to make measurements of or then clearly the above density matrix cannot be 
distinguished from the maximally mixed state 

_ 1 oA 
''"““2(01)’ 



which has entropy In 2. This is precisely the same entropy we would attribute to the two spin 
system had we taken the second spin to be a true extra degree of freedom, although it is clearly 
not the standard meaning of entanglement entropy. But had we started with two real spin degrees 
of freedom, we can also think of the entropy resulting from tracing out the second spin as the 
ignorance introduced when we stop making measurements of operators of the form ® (^lyzi 
where G {I 2 , That is, we are only allowed to make measurement of the first spin using 

operators O G {o'j. ^ ^ (8) 12 }- This is of course the familiar notion of keeping only local knowledge 
within a chosen “region”. What we attempt to illustrate with this trivial example are two points: 
first, that entanglement can be viewed as ignorance after dumping some observables from being 
measured; second, which trivially follows from the hrst assertion, is that the choice of measurable 
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observables directly affects the final value of the entropy. This is precisely the point made in Ref. 2 
over the effect of the choice of the operator algebra Ay on the value of the entanglement. The 
entanglement entropy following from the electric center choice is the ignorance introduced from 
dumping the loop magnetic operators at the boundary of the region. This is the direct analogue 
of the above entropy for the single spin system where we dump from our measured set. To 
reiterate, the procedure of deleting off-diagonal components in the density matrix after expressing 
it in the eigenbasis of the electric link operators Lf is indeed simply equivalent to dumping magnetic 
operators non-local with the boundary Lf. Whether this is considered classical contribution to 
the entropy, as asserted by Casini is perhaps more or less a matter of taste. 

As has been emphasized in Ref.2 therefore, an ambiguity, or more precisely, some artihcial 
choice one has to make in dehning the entanglement entropy is not only special to gauge theories, 
but to all kinds of theories, when we make different choices of sets of operators that we perform 
measurements. The major difference of a gauge theory compared to a local scalar held theory 
however, is that there exists, in the case of say the local scalar held theory, a natural or canonical 
choice of the operators that we keep, which dehnes a notion of spatial entanglement, and that such 
a natural choice is absent in a gauge theory. 

The crucial question, therefore, boils down not so much to which is a correct choice of operator 
algebra Ay, but rather, which is a useful choice that would give us useful information about the 
system under scrutiny. 

Here, we recall that entanglement entropy (or more precisely, mutual information) is a measure 
of correlations between diherent degrees of freedom. Of course, strictly speaking the correlation 
following from (13) is trivial - states not satisfying the condition simply does not exist. Never¬ 
theless, realistically, in an experimental setting, it is not the Hilbert space that we are directly 
probing, but the expectation values of (gauge invariant) observables. Therefore, from the perp- 
sective of an experiementalist, which do not have a clue of the redundancy of the Hilbert space, 
would be fascinated by the agreement between the results of al and a^. Therefore it is perhaps 
too immature to discard the idea that the entanglement entropy of the electric center is just as 
real an entanglement entropy as it could be. 

We note also that the dramatic difference between the electric centre and magnetic centre simply 
follows from the fact that the ground state behaves as a direct product state if we use a basis that 
consists of the eigenstates of the magnetic flux loop operator. We note that this is precisely the 
same basis for the Z 2 symmetry protected topological (SPT) phase, recalling the “duality” relation 
between SPT and gauge theory as hrst described in Ref.5 and elucidated in detail in Ref.6. It is 
therefore not surprising that the topological entanglement entropy should vanish. 


IV. GENERALIZATION TO NON-ABELIAN THEORIES 


In non-Abelian gauge theories, the degrees of freedom again live on edges, and are labeled by 
group elements of a gauge group G. One needs a convenient representation of the gauge invariant 
Hilbert space in order to discuss the splitting of space. This is well-known in tensor network theory. 
The idea is that one can perform a discrete generalized Fourier transform on the states on each 
edge, as follows^ 

|/i,a,6) = (16) 

g V I I 

whereas the inverse transform is given by 

\ 9 ) = (17) 

Note that we abused the notation /i for both an irreducible representation of G and the dimension 
of the representation. 


Under a gauge transformation therefore, the new basis transforms as 

9-^1 I h) in in ^ 9-^1 1/^; ^)in {,9 )ac |/^) C, 6) in 

9^1 1 99) out I 9i9 )out ^ \pi Q) 99) out \pi (^1 C )out P^{,9')cbi (1^) 

where repeated indices are summed. 

The ground state can now acquire a gauge invariant representation as we translate the ground 
state into this representation. Recall that the ground state of the theory is given by equation 
(5). While we have emphasized that configurations related by gauge transformation in a gauge 
theory is strictly identihed, one could think of (5) as a non-gauge-fixed representation, to make 
the ground state appear explicitly gauge invariant. 


Now, if we just follow through the change of basis, the ground state would look like the following 

(19) 


i;Evertices g /Glinks 


appropriate orientation ) 


It’s not hard to understand what the projectors do. To make it explicit, we focus, without loss 
of generality, on a particular 3-point vertex, and let us assume that all links are pointing outwards. 
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Then the action of the projector gives 


E E 




hlh2h3 


l/il, ai, Cl) |/i 2 , 02 , C 2 ) |/i 3 , 03 , {9)c2b2p^^ (^)c 363 


where 


J2p''"^3)c^b^P^%9)c2b2p^^{9)c^b^ = ^(o;hTEr)'a; 


|G|, 

C1C2C3 / ^616263 ’ 


h3 


( 20 ) 

( 21 ) 


where the ool^"'} are the intertwiners that map prodnct representations of ui 0 U 2 to representation 
The factor of l/ps follows from the identity 


1 

M 


y^P^{9)abP^{9~%cl = -SadSbc. 

„ P 


( 22 ) 


The statement is simply that by snmming over g, the snrviving terms correspond to the singlet 
representation by decomposing the prodnct representation of Hi® 112 ® hs, simplicity we 

assnme there is only one fnsion channel to the singlet. Otherwise we will have mnltiple versions of 
cj, each corresponding to a particnlar fnsion channel. Had the vertex been endowed with a valency 
greater than 3, the end prodnct wonld be a series of fnsions that take the representation to the 
singlet representation. Therefore, we have 




E 

fJ-X 


/il/X 2 

1^ 


^clc 2 lfy^bXbT 1/^2, 02 , C 2 ) \pi, 02 , C 2 )\P3, 03 , C 3 ) ' ' ' , (23) 


and if the system consists of this single vertex, we shall have o* = bi that are also snmmed over bnt 
will not have any term in the neglected part ’•••’. We thns recover the well-known fact in tensor- 
network theories, that the gange invariant Hilbert space consists of all conhgnrations specihed by a 
set of {/ij} at each links, snch that when mnltiple links meet at a vertex, there mnst exist a fnsion 
channel snch that they fnse to the trivial representation. If more than one snch fnsion channel 
exists, we wonld have to specify the particular channel. Therefore a state is specihed by 


l^) = IW,{7r}), (24) 

where we take to denote the particular fusion channel to the singlet at a given vertex Uj. The 
internal indices are not gauge invariant; hence, they are not needed to specify the Hilbert space. 

Now the question is, if we were to cut the system into two pieces in some sense, how would 
we attribute entanglement entropy to these states? On the one hand, the internal indices are not 
gauge invariant and thus unphysical. On the other hand, we have no other way to specify the 
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states without referring to them, and the particular contraction of the link internal indices across 
vertices again make it ambiguous to dehne a product space, and it appears that we are no better off 
than the previous basis. At this point, we can draw experience from our discussion in the previous 
section and also the case of Abelian theories. Recall in the case of Abelian theories, the outcome of 
measurements of “spins” connected by a vertex are correlated due to the gauge invariant constraint. 
While we argued that such a entanglement might be trivial, it is still of interest to quantify these 
observed correlation, by discarding a set of observables from the set of observables that would be 
measured, and ask for the density matrix with the largest amount of entropy that would agree 
with the original density matrix over the set of measured observables. Phrased in these terms, it is 
clear what measurements we would like to keep making. We would like to be able to measure /Xj, 
which are gauge invariant variables, and we would hnd uncanny correlation in these measurements 
of these representation labels /Xj when the links meet at a vertex. As in the case of picking electric 
center of Abelian theories, any such “entanglement” entropy has to arise at the boundary, where 
we give up measuring observables not commuting with the measurement of /i^. In non-Abelian 
theories, this would again boil down to Wilson loops, much like what happens in Abelian theories. 
To compute the resultant entropy that characterizes our ignorance, we would like to start with 
computing the probability density of obtaining a particular set of outcomes {/i} for links taken to 
be in region A. In the interior region, these outcomes have no dependence on dofs in A, and only 
at the boundary they are controlled by links in A via Gauss’s law. Once we obtain the probability 
density of individual outcomes for these boundary links (picture exactly like the electric center 
in^), we would then be able to recover the entropy by 


s = - lnp{/ib}i, ( 25 ) 

i 

and this would be understood as the entanglement entropy. Again let us emphasize this is precisely 
the same procedure as in Ref.2 of obtaining a density matrix prepared in a basis that diagonalizes 
the electric observables, and then “deleting” off diagonal entries before computing a “classical” 
entropy arising from the distinct blocks in the resulting mixed density matrix. 

Let us begin with a simpler example, in which the whole system consists of again precisely one 
vertex connecting three links. We will take one link, labeled by /is, as belonging to region A, and 
the other two to A. The entanglement of integrating out A can be computed using the form of the 
ground state as given in (23). The reduced density matrix is 


PA 


Mi 


P1P2 


I //1//2M3 |2| /ii/i2M3 |2 

\^aia2a3 \ l^ciC2C3 I 


\p3, 03) C3)(/i3, 03 , C 3 I -f • • ■ 


(26) 


where • ■ • denote components not diagonal in /is. Therefore, the probability density of measuring 
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/i 3 is given by 


P{P-'3) 


MI1M2, 

Atl®A‘23A‘3 




Ml)M2 




|GP 


where in the hrst equality above we have made use of the fact that 


(27) 


Y,\<SS\" = ^‘3■ ( 28 ) 


Also, the fusion product of of /ii®/i 2 is supposed to contain ps, but by inserting corresponding 
to the fusion matrix, the sum is effectively unconstrained, and hence the second equality. But keep 
in mind, since we assumed the coupling of representations in this case is multiplicity free, here in 
fact which is either 0 or 1. 

Now using the fact that for hnite groups. 
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Y,ax“(9) = \G\s,.„ 


(29) 


where x'^id) characters of the representation /i. Substituting into (27), we end up with 

P(A3) = |§|, (30) 

which is an answer that makes a lot of sense. To start with, note that = 1- This answer 

suggests that despite not being to measure in detail the internal states of each given representation 
p, the probabilities are in fact sensitive to the dimension of the representations. 

Now having done the warmup exercise of a single vertex, we will come to the meat of our 
discussion. We will consider some arbitrary region in a system with many vertices and links, and 
try to dehne an entanglement entropy with the regions beyond. 

Consider region A surrounded by a set of links forming the boundary. The links touching the 
boundary belonging to A would depend on the state that the boundary links of A take due to 
the Gauss constraint, and we would again like to dehne an entanglement entropy along the same 
vein as we have done for the single vertex case. Here, instead of integrating out everything in A, 
we could equally ask for the probability of hnding the set of protruding links in A touching the 
boundary of A in a particular conhguration {/iz}, and that probability would be used for dehning 
the entanglement entropy. The computation pretty much parallels the single vertex case above, 
with the only new ingredient being that each leg lying along the boundary of A has to connect 
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to two different vertices, and that they form a closed loop, such that one can expect some extra 
global relation between different touching links in A, again in complete analogy with the Abelian 
case. 

Let us assume that there are L boundary links forming a loop, with orientation such that 
altogether it is an counter-clockwise loop. The L boundary links are separated by L touching links 
in A. We note that the external links and the boundary links are meeting at trivalent vertices for 
simplicity and concreteness. Let us allow the links to point toward the loop. Now we would like 



FIG. 3. 


to compute the probability of finding the set of L external links in any particular conhguration 
denoted by {i^j\j = 1, • • ■ L}. The representations residing on the boundary of A would be denoted 
as {fii\i = 1, • • ■ L}. 


Let us hrst write down the state 

L 




ly-i 


i 1 / C'iCl-i 




1^1 r/9 TT A / \Q\P^A9i-l9i )cibi\f^i,bi,Ci) ,^,pA9j)mjnjWj,'^j,nj), (31) 

{gi} i ^ ' j ^ ' 


where the indices i,j are taken to be cyclic modulo L, i.e., i = L = 0 (mod L). 

By summing over everything in A, we can construct the “reduced” density matrix as follows, 

{gi},{hi},{mj,nj} i 

^ {9j)mjnjP ^j) (z/j, rUj, nj \ + • ■ • (32) 

i 

where we only write down the diagonal components explicitly, and all off-diagonal components are 
implicitly included in • • •. 

Now we would like to compute the probability of finding the external legs in a particular 
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configuration {r'i}. The probability would involve the sums 


En 

{/^*} ^ 






-1 


hih -\)= 


9i-l9i 


(33) 


This product of delta functions imply that 

^ ^ Qik ^ hi, Vi', 


(34) 


for any k E G. Hence, summing over {gi} and {hi} independently can be traded for summing over 
{qi} and k. We are then left with 


Pi{G}) 


fl-'-VL 


k,{9i} * 


Ui-'-VL 

\G\^ 


E fJXl ATX2 

1/11^2 a:ii/3 

k,{xi} 


. .. V' 

XL-2^L^ 


XL- 


^{k) 


Ui-'-Vl 


E 


prxi pTX2 . . . pTXL-l=l 
-'’1/11/2 XIZ/3 XL-2h'L 


{xi,...XL-2} 


(35) 


The arises from p'^'{9i)miniP^^*{9ik)mini- These probabilities clearly sum to 1, again using 



(36) 


In the above, we are considering a region A that essentially does not contain an interior, and all 
its links he at the boundary. To consider a more general region H, there would generally be more 
links at the boundary compared to the v links from A touching the boundary of A. Nonetheless, 
as one can see from the above, as we sum over the /r's, eventually, all that imposes is 

97\9i = K\hi, (37) 

despite that now the index i would run from 1 to A, where A > L in general. Therefore, when it 
comes to the probability distribution of hnding the external u links in a particular conhguration 
P{{ui}), the effect is entirely identical to the above where the interior is trivial and A = L. 

Let us also remark that if we compare these results with Ref.8, there is a reduction in entropy, 
precisely due to the fact that the internal states are not gauge invariant and are not detectable 
even experimentally. 


14 





V. CONCLUSION 


In this note, we attempt to clarify what the purported ’’electric entanglement entropy” is 
measuring, and why it is still an interesting quantity. In the Kitaev model where the physical 
Hilbert space contains non-Gauss projected states, one recovers exact electric magnetic duality i.e. 
there is no difference between the so called magnetic center and electric center. We generalize our 
considerations to non-Abelian gauge theories. 

During the preparation of the manuscript, we become aware of the work,® which partially 
overlaps with our current discussion. In particular, the discussion in® emphasizes the fact that 
the gauge invariant physical spectrum is conveniently embedded in a larger Hilbert space, which 
is given by the value taken by each link without worrying about gauge hxing/Gauss’s constraint. 
That is precisely the Hilbert space of the Kitaev model that we explored at the beginning of the 
note. 
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